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Abstract. We report on a search for low-energy neutrino (antineutrino) bursts in correlation 
with the 8 time coincident events observed by the gravitational waves detectors EXPLORER 
and NAUTILUS (GWD) during the year 2001. 

The search, conducted with the LVD detector (INFN Gran Sasso National Laboratory, 
Italy), has considered several neutrino reactions, corresponding to different neutrino species, 
and a wide range of time intervals around the (GWD) observed events. No evidence for 
statistically significant correlated signals in LVD has been found. 

Assuming two different origins for neutrino emission, the cooling of a neutron star from 
a core-collapse supernova or from coalescing neutron stars and the accretion of shocked 
matter, and taking into account neutrino oscillations, we derive limits to the total energy 
emitted in neutrinos and to the amount of accreting mass, respectively. 
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1. Introduction 

The analysis of the data collected in coincidence by the gravitational wave bar detectors 
EXPLORER and NAUTILUS during the year 2001 (Astone et al., 2002) shows an excess (8 
events against 2.6 expected from the background) when the two detectors are favorably oriented 
with respect to the Galactic Disc. Moreover, this result comes from the present day most sensi- 
tive experiments for the detection of gravitational wave bursts and a search for neutrino bursts in 
correlation with the 8 GWD events is, therefore, appropriate. 

A few astrophysical transient sources are indeed expected to produce associated bursts of 
neutrinos and gravitational waves. It is well known that most of the energy (99%) released in 
the gravitational core collapse of a massive star is carried away by neutrino originated both from 
the matter accretion in the shock and from the cooling of the proto-neutron star (see for example 
Burrows et al., 1992). Depending on the collapse dynamics, some fraction of the total energy is 
emitted in GW ( Thorne 1988, Muller 1997 ), asymmetric supernovae in our Galaxy being the 
best candidate sources for GW bar detectors. Two coalescing neutron stars would also constitute 
a source for both neutrinos and gravitational waves. From the point of view of GW emission, it 
is likely that the merging event would produce powerful gravitational wave bursts, and, even if 
the physics of the merger is not known, there are estimates that, for binary systems of large mass, 
coalescence waves are likely to be stronger than the inspiral ones. Some amount of the kinetic 
energy is converted in thermal energy so that the hot remnant would probably emits thermal 
neutrinos. 

The search for a neutrino burst associated to the events detected by the GWD EXPLORER and 
NAUTILUS in 2001 has been performed with the LVD apparatus, operating in the INFN Gran 
Sasso National Laboratory (Italy) since 1992 with the main purpose of searching for neutrinos 
from gravitational stellar collapses within the whole Galaxy. 

The paper is planned as follows: in Sect.2 we briefly describe the LVD detector, and we 
explain the selection of the data. In Sect. 3 we present the results of the analysis: a time inter- 
val spanning from 12 h preceding each of the 8 events up to 12 hours later, has been scanned, 
searching for any excess over the statistical fluctuation of the background. Further, a search for a 
v signal coincident in time with every event has been performed. We conclude in Sect.4, where 
we discuss the results of the search, taking into account v oscillations, and considering the two 
processes for v emission, i.e., cooling and accretion. Since we do not find any neutrino burst 
candidate associated with the 8 GWD mentioned events, in the following assumed scenarios we 
derive upper limits: 

- on the neutrino flux, without reference to any particular source, 

- on the total amount of energy emitted in neutrinos, in the cooling case, 
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2. The LVD Experiment and the Data Selection 

The Large Volume Detector (LVD) in the INFN Gran Sasso National Laboratory, Italy, consists 
of an array of 840 scintillator counters, 1.5 m 3 each, interleaved by streamer tubes, arranged in a 
compact and modular geometry (see Aglietta et al., 1992, for a detailed description). The active 
scintillator mass is M = 1000 t. The experiment has been taking data, under different larger 
configurations, since 1992, and reached its final one during 2001. 

There are two subsets of counters: the external ones (43%), operated at energy threshold £/, 7 
MeV, and inner ones (57%), better shielded from rock radioactivity and operated at £/, ^4 
MeV To tag the delayed y pulse due to n-capture, all counters are equipped with an additional 
discrimination channel, set at a lower threshold, £/ - 1 MeV. 

The main purpose of the telescope is the detection of neutrinos from gravitational stellar collapses 
in the Galaxy. In the following we will focus on v reactions with free protons and 12 C nuclei, 
constituting the bulk of the expected signal and having the best signature in the detector, namely: 

- (1) inverse y6-decay: v e p, e + n, observed through a prompt signal from e + above threshold £/, 
(detectable energy Ea - E Yg - 1.8 MeV +2m e c 2 ), followed by the signal from the np,dy 
capture (E y = 2.2 MeV), above £/ and with a mean delay At 180 fis. 

- (2) Vi and v t neutral current interactions with 12 C: v ( - 12 C, V,- 12 C* (z = e,fx, r), whose signature 
is the monochromatic photon from carbon de-excitation (E y =15.1 MeV), above & h . 

- (3') v e charged current interactions with 12 C: v e 12 C, 12 N e~, observed through two signals: 
the prompt one due to the e~ above £/, (detectable energy Ed E Ve - 17.3 MeV) followed by 
the signal, above £/,, from the fi + decay of 12 N (mean life time t = 15.9 ms). 

- (3") v e charged current interactions with 12 C: v e 12 C, 12 B e + , observed through two signals: 
the prompt one due to the e + (detectable energy Ed Ey e - 14.4 MeV + 2m e c 2 ) followed by 
the signal from the /T decay of 12 B (mean life time t = 29.4 ms). As for reaction (3'), the 
second signal is detected above the threshold £/,. 

After being subjected to a preliminary process to reject muons, the raw data are grouped in three 
different classes, with specific signatures to tag the different described reactions: 

1. IBD class (inverse beta decay): pulses with E4 > £/, followed by a delayed (At < 1ms) 
low energy (E& > £;) pulse in the same counter. The efficiency in tagging the n-capture is 
e = 60% for the core counters, e = 50% for the whole detector; 

2. NC class (neutral current): pulses with 11 MeV < E d < 17.5 MeV, the efficiency in tagging 
the y from 12 C de-excitation being e = 55% 

3. CC class (charged current): two pulses with E d > & h MeV, within At < 150ms, in the 
same counter. The efficiency in tagging the 12 N and 12 B decay are e = 85% and e = 70%, 
respectively. 
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3. The Analysis 



The LVD detector is sensitive to neutrino bursts from core collapse supernovae within the whole 
Galaxy (Aglietta et al., 1992). The scintillator counting rate is continuously monitored: all the 
events are examined on-line on the basis of their time sequence. Neutrino burst candidates are 
identified as clusters of scintillator counter pulses with an imitation frequency less than a prede- 
fined threshold (Fulgione, Mengotti and Panaro, 1996). During the year 2001, no neutrino burst 
candidate has been evidenced, thus allowing to conclude that no v signal from gravitational stel- 
lar collapse in the Galaxy has been detected (Aglietta et al., 2003). 

However, the absence of candidates in the LVD detector taken alone does not preclude the pos- 
sibility of positive effects, when combining it with another detector, since the joint measurement 
allows to increase the sensitivity. The analysis in correlation with the 8 candidate events has then 
been conducted, in four steps described in the following. 



3.1 . Step 1. Check of the detector stability. 



First of all, the LVD detector performance at the occurrence of the 8 GW events (see the list in 
Astone et al., 2002) has been checked by studying the behavior of the counting rate in a 24 hours 
interval around the time of each of them. 
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Fig. 1. Distributions of the deviations of «is (see text), with respect to («is), for events of IBD 
class (top) and NC class (middle), and nis distribution, for CC class (bottom). The superimposed 
dotted lines are the fits to zero mean and unit-width gaussians (IBD and NC) and poissonian 
(CC). 



For each event, the 24 hr average of the number of counts every 15 minutes, («is), is evalu- 
ated: including all the 8 events, we show in figure^the distributions of the deviations of n^, with 
respect to (nis), for events of IBD class (top) and NC class (middle); for CC class (bottom), due 
to the small < nis >, we directly show the nis distribution. They have been fitted with zero mean 
and unit-width gaussians (IBD and NC) and poissonian (CC): they are shown superimposed in 
the same figures, together with the resulting reduced^ 2 values. The LVD counting rate, for all 
the 8 events and all the data classes, is then well understood in terms of poissonian statistics: this 
sets a firm base for the following steps. 

3.2. Step 2. Search in a sliding window. 

The search for a possible v burst has been conducted in a 24 hours interval T around the occur- 
rence of each of the 8 events. The 8 intervals have been scanned through a "sliding window" of 
variable duration: more in detail, they have been divided into Nst = 2 ■ £ — 1 intervals of different 
duration 5t, each one starting at the middle of the previous one. The multiplicity distributions of 
clusters (i.e., the number of events within each 6t) have then been studied for the three classes 
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of data and for 5t = 1,5,10, 20, 50, 100 s, and have been compared with the expectations from 
poissonian fluctuations of the background. 




Fig. 2. Distributions of cluster multiplicities, for IBD class events (top), NC class (medium), CC 
class (bottom), together with expectations from poissonian fluctuations, in the case of GW event 
n.5. 



We show as an example the case of GW event n.5: the distributions of cluster multiplicity, 
for events of IBD class (top), NC class (middle) and CC class (bottom), in the case of the 6 dif- 
ferent 5t, can be seen in fig.|2] together with the expectations from poissonian fluctuations of the 
background, the relative reduced^- 2 values ranging from 0. 1 to 1 .2. The agreement between data 
and expectations holds also in the case of the other seven events. This, together with the check of 
the poissonian probabilities associated to each measured multiplicity, in each class and for each 
event, allows to state that there is no evidence for any detectable v signal in correspondence of 
any of the considered events. 

3.3. Step 3. Search in a fixed window. 

The search for a v signal in coincidence with every GW event has been further conducted us- 
ing a "fixed window" centered at the time of each of them. In particular, for each data class, 
we compare the number of pulses (iVd), recorded during time windows of different duration St, 
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centered on each event time, with the average number of pulses expected from background, N^k- 
The value of Nbk is evaluated by using the rate in the 24 hours around each event, excluding the 
contribution of the central portion of time to avoid the contamination due to a possible signal. 
Results corresponding to St — 1,5, 10, 20, 50, 100 s, for each of the 8 events, are summarized in 
table^EJfor IBD and NC classes, respectively. The table relative to CC class is not shown since 
in all the cases Nd = 0, and Nbk < 3 ■ 10~ 2 . 



Table 1. IBD class: number of events detected in coincidence with the 8 GW events, for different 
durations of the time window (St), compared with the expectations from the background. The 



effective LVD mass, M, at the time of each event is also shown. 
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Table 2. NC class: number of events detected in coincidence with the 8 GW events, for different 
durations of the time window (St), compared with the expectations from the background. The 
effective LVD mass, M, at the time of each event is also shown. 
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The differences between Nd and N\,k are within the statistical fluctuations, for all data classes 
and for all the events. The most significant effect is observed in correspondence of GW event 
n.5: 11 pulses detected against 3.4 expected, when using IBD class data and 5t = 100 s. Taking 
into account the number of trials (8 x 3 x 6), the associated chance probability is P — 0.03. In 
order to check the consistency of such an effect with a physical one, we complete the coincidence 
analysis with the study of the time distribution of both high and low energy signals. 

3.4. Step 4. Time distribution of pulses. 

We have studied the time distribution of LVD pulses around each GW event. Figure[3]shows such 
a distribution for each of the 8 events (and for their sum), for IBD data (full line) and NC data 
(dash-dotted line), (r = corresponds to the time of the GW event): as for all the other classes, 
no particular time structure is present. 
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Fig. 3. Time distribution (bin=5s) of the detected pulses (full line for IBD class data and dash- 
dotted line for NC class data) around the corresponding GW event time (f = 0); the bottom panel 
shows the distribution for the 8 events taken together. 



Finally, for IBD class data, we have also checked the time distribution of secondary pulses 
(i.e., those possibly due to neutron capture) with respect to the prompt ones. The measured dis- 
tribution is uniform and compatible with the one expected in the case of pure background, where 
delayed and prompt signals are uncorrected and the distribution of the differences in time is 
flat (On the contrary, if the pulses were due to v e interactions with protons, the distribution of 
time delays should show an exponential behavior, with t ~ 180/is, corresponding to the average 
capture time of neutrons in the LVD counters). 

4. Calculation of Upper limits 

4.1. Upper limits on neutrino fluence 

No evidence for any statistically relevant signal in LVD, in the three considered reaction channels 
(corresponding to different neutrino species) and over a wide range of time durations, has been 
found in correspondence of any of the 8 excess events detected in coincidence by NAUTILUS 
and EXPLORER. 

In the absence of any v signal, we calculate 90% C.L. neutrino fluence upper limits at the de- 
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tector without assuming particular energy spectra, i.e., on mono-energetic neutrinos at different 
energies: 

^= M.N t .aiE v yeiE v) (1) 
where: A^o is the 90% c.l. upper limit on the number of LVD signals per GW event, obtained 
following Montanet et al. 1994 in the case of Poisson processes with background. The considered 
background value represents the total number of expected background events for all the 8 GWD 
events; the signal value, as well, is the total number of detected events for all the 8 GWD events. 
M is the detector active mass in ton (summed over the 8 events); N t is the number of targets per 
ton (either protons or 12 C nuclei); e is the detection efficiency; cr(E v ) is the appropriate cross 
section. Results are shown in table[3]for St = 20s and 100s. 



Table 3. Fluence upper limits (90% C.L.) for neutrinos of different energies, obtained from IBD 
and NC classes of events. 
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This model independent fluence can be used to test a specific model (characterized by 
Q> test (E v ) ) by performing the convolution: 

MOOIteV 

x= ———dE v (2) 

J\0MeV ^tesA^v) 

and if the quantity x is less (greater) than 1.0, the model predicts more (fewer) events than the 
event limit N90, and is therefore excluded (not excluded) at the 90% confidence level. 

4.2. A model dependent interpretation 

We can discuss the result of the search within two possible simplified scenarios for neutrino 
production, namely (i) thermal emission (which we will call "cooling") and (ii) non thermal 
emission (which we will address as "accretion"). 
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We assume that the 8 events are due to a unique kind of source and that the distance to the 
hypothetical sources is d — 10 kpc, since the 8 GWD events are consistent with a galactic origin. 
Concerning neutrino oscillations (Dighe and Smirnov, 2000; Takahashy et al., 2001; Aglietta et 
al., 2003), due to the unknown v oscillation parameters, i.e., |t/ f 3| 2 and the mass hierarchy, we 
consider four different cases: 

1 - normal mass hierarchy with adiabatic transition at the high density (H) resonance {NH a d)\ 

2 - normal mass hierarchy with non-adiabatic transition at the H resonance {NH non - a d); 

3 - inverted mass hierarchy with adiabatic transition at the H resonance (7// fl£ /); 

4 - inverted mass hierarchy with non-adiabatic transition at the H resonance (IH, m „- a( i). 

More details on the effect of neutrino oscillations in the supernova neutrino signal in LVD can 
be found in (LVD coll., 2004). 

We refer to Appendix A for details on the calculation of the number of events in each detec- 
tion channel. 



4.2.1. Cooling 



In the simplified "cooling" process we are considering (for example in the case of a newly formed 
neutron star), neutrinos of every flavor are emitted by electron-positron annihilation (e~e + , v,v,, 
i = e, fi, t) with thermal spectra, that is, we are assuming zero pinching. We assume exact equipar- 
tition of the total emitted energy E B among all neutrino flavors (E Yj = f Vj E B with f v . = 1/6) and 
a hierarchy of the mean temperatures of the different flavors T Yx > Ty e = T Ve (x = //, t). The 
characteristics of this emission process and the numerical values used are summarized in tab. |4] 
and in the Appendices. 



Table 4. Characteristics of the considered processes of v emission and numerical values used in 
the limits calculation. 
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Fig. 4. 90% c.l. limits to the total energy emitted in neutrinos, in the cooling process, during 
20 s, assuming pure Fermi-Dirac spectrum. The dashed line represents the IH al j case, while the 
other three cases are represented by the solid curve. The limit obtained through NC events is 
represented by the dash-dotted line. 

The limits obtained using IBD events - which are the most stringent (see the appendices for 
the calculation) - in terms of total emitted energy Eg, are shown in fig.|4]as a function of 7V e , in 
the case of exact energy equipartition among all flavors, with T Vx - 1.57V and 5t = 20s. 1 The 
dashed line represents the IH a d case, while the other three cases are not distinguishable among 
them and are represented by the solid curve. For the sake of completeness, we shown in the same 
figure the limit obtained through NC events (dash-dotted line) which is independent on v flavor. 
As an example for the case of cooling, we show in the same figure the total energy expected to 
be emitted in two of the most probable sources: a new-born neutron star (Keil, Raffelt and Janka, 
2002, and references therein) empty box, and colliding neutron stars, full box (Ruffert and Janka, 
1998). 

4.2.2. Accretion 

During matter accretion in a neutron-rich ambient, another process appear to generate v emission 
(see e.g. T.J.Loredo and D.Q.Lamb, 2002): e ± pairs, created in the accreting matter, are captured 
by protons and neutrons and the resulting v fluxes have the following characteristics: (i) only v e 

1 Using the values: f Ve = 1/5 and T Vx = 1.37? , the limits weakens of a factor < 2. 
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and v e are emitted; (ii) the thermal energy spectrum is multiplied for a E 2 V term accounting for 
the capture cross section energy dependence. The characteristics of this emission process and the 
numerical values are shown in tab.0] The obtained limits (see the appendices for the calculation) 
in terms of accretion mass multiplied by the neutron fraction of accreting matter, Y„, are shown in 
Figure|5J as a function of 7\> e . The solid line represents the limits obtained using IBD interactions 
in the case of NH and lH non - a d- Since for IH ac i only a minimum part (\U e 3\ 2 ) of the original v e 
flux will interact as v e , we do not use the IBD events. It still remains valid the one obtained with 
NC events, represented by the dash-dotted line. 




Tv e [MeV] 



Fig. 5. 90% c.l. limits to the total accretion mass, multiplied by the factor Y n , as a function of 
Ty e . The solid line represents the limits obtained using IBD interactions in the case of NH and 
lH non -ad, the dash-dotted one refer to limit derived through NC events. 



As an example, we show in the same figure (full box) the amount of accreted matter (multi- 
plied by Y„ = 0.5) expected for coalescing neutron stars (Ruffert and Janka, 2001). 



5. Conclusions 

We have conducted a search for low-energy antineutrino (neutrino) bursts with the LVD detector 
in coincidence with the 8 event excess found by the gravitational waves detectors EXPLORER 
and NAUTILUS during the year 2001. 

We have found no evidence for any statistically relevant signal in LVD, in three different 
reaction channels (inverse beta decay, charged current and neutral current interactions with 12 C) 
corresponding to different neutrino species, over a wide range of time durations, for any of the 
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8 events. Consequently, we have derived 90% fluence upper limits to antineutrino and neutrino 
emission from an average GW event, at different energies in the range of sensitivity of the LVD 
detector. 

We have then related the result of the search with two possible simplified models for neutrino 
emission, i.e., "cooling" and "accretion", deriving limits, on the one side, to the total energy 
emitted in neutrinos at the source, and, on the other, to the amount of accreting mass. Assuming 
a source distance d — 10 kpc, possible candidates as new-born and colliding neutron stars have 
been excluded by this analysis. This makes even more challenging and interesting the search for 
a likely astrophysical source for the reported GWD events. 
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Appendix A. Limits calculation: inverse p decay interactions. 

The number of v e interactions due to inverse beta decay in a detector is given by: 

r°° dN- 

N' e B v D = M-N p - —^■<r(E-^-s n -s(E d ,E th )dE % (3) 
Jq uhy e 

where: M is the detector active mass in ton; N p = 9.36 • 10 28 is the number of free protons 
per ton; e„ is the neutron detection efficiency; e(Ea, E t / t ) is the e + detection efficiency; Ed = 
Ey e - Q + 2m e c 2 , with Q = M n + m e - M p = 1.8 MeV is the positron detectable energy; E,h is 
the detector energy threshold; cr(E Ye ) is the cross section (Vogel and Beacom, 1999; Strumia and 
Vissani, 2003); is the antineutrino energy spectrum. 
In the case of cooling process it is: 

dE- Ve And 2 In 4 v " ^ ' 

where d = lOKpc is the assumed source distance; 

Eg = Yji fv.Es is the total energy emitted in neutrinos; 

F v is the term accounting for different v oscillation scenarios 2 : 



- for NHaH\ NHnon-ad', IHnon-aj: 
fve it, i2 E 2 /> 

r 4 1 ^ 11 i+e E/T> e + r 4'— l +e E/T, 

V e V x 

2 NC data are not affected by oscillations. However, the limits from IBD data stay almost the same 
even assuming that MSW oscillations are completely absent. In fact, on accounting for vacuum oscillations 
we get Pg g = 1 - sin 2 20 i2 /2 ~ 0.6 in all scenarios; this is practically the same value implied by MSW 
oscillations in the scenarios NH and IH„„„ ad _, P- ee = cos 2 6[ 2 ~ 0.7. 
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- for IH ad : 

with: \Uei\ 2 ~ cos 2 6 n = 0.67, |[/ e2 | 2 ~ s/n 2 #i2 = 0.33, \U e3 \ 2 > 10~ 4 for the adiabatic case and 
\Ue3\ 2 < 10~ 6 for the non adiabatic one (M. Apollonio et al., 1999) (J.Bahcall and C.Pena-Garay, 
2003). 



In the case of mass accretion process it is (Loredo and Lamb, 2002): 
dN v 1 



A a ■ Y n ■ M hot ■ F' (7) 



dE % And 2 
where: 

Mhot is the mass of hot emitting material; 
Y n is the neutron fraction; 

A a = l+3 / A ; r ° c ,,, tI^t, with g A = 1.254, o-o = 1.7 • 10 44 cm 2 ; 
and with respect to i 7 ': 



F; = |t/ el | 2 - — — (8) 

- for IHad: 

F' 9 = |t/ e3 | 2 £ ' - (9) 

Appendix B. Limits calculation: neutral current interactions 

The number of interactions in the detector due to the neutral current is given by: 

N r r°° dNy. dN v 

N^ c =M-N c -e c - ["TF 1 ^) + -nr<r(E Vi )]dE (10) 

J15.1 MeV dt,y i atL Vj 

where: Nc = 4.24 • 10 28 is the number of 12 C nuclei per ton; sc is the detector efficiency for 15.1 
MeV gamma; a-(E v ) is the cross section (M. Fukugita et al., 1988). 
The neutrino energy spectrum in the case of cooling process is: 

dl% = d]% = _Eb_ 120 

dE- Vi ~ dE Vi ~ And 2 ' ln A ' '' 

with Fi = ji- 1+ ^ /r ,. f + 2 • |r 1+e £/r v . v 

while, for the mass accretion case, we considered all the events as v e and we used: 

dNy. dN v 1 E 4 

a^r^tr^ Aa - Yn - Mhor T^z (12) 
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